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highest  energy  or  intensity;  the  spectrum  of  the  stiff  soil  sites  is  slightly 
less  than  for  the  rock  sites;  and  the  spectra  of  the  soft  clay  and  sand  sites 
and  the  deep  cohesionless  soil  sites  are  almost  the  same  and  are  lower  than 
those  for  stiff  soils.  However,  in  the  low-frequency  range  of  0.006  to  2.5  Hz, 
the  reverse  exists:  the  soft  sites  indicate  the  highest  energy,  the  deep  co¬ 
hesionless  soil  sites  are  next,  the  stiff  soil  sites  are  third,  and  finally, 
the  rock  sites. 

A  qualitative  comparison  was  made  of  the  spectral  shapes  of  the  PSD 
calculated  in  this  study  with  the  acceleration  response  spectra  (ARS)  of 
Kiremidjian  and  Shah  (1978),  Seed  and  Idriss  (1971),  and  Seed,  Ugas,  and  Lysmer 
(1976).  There  is  general  agreement  in  the  shapes  of  both  types  of  spectra, 
except  the  relative  amplitudes  for  the  rock  sites  for  the  ARS  were  less  than 
those  of  the  PSD  in  the  high-frequency  range.  The  difference  in  relative  ampli¬ 
tudes  for  the  rock  PSD  and  ARS  is  due  principally  to  the  lesser  number  of 
records  used  in  the  data  analysis  for  the  latter.  The  number  of  records  for  the 
rock  sites  was  28  for  the  ARS  and  56  for  the  PSD.  Another  difference  is  that 
the  spectral  shapes  of  PSD  functions  show  more  peaks  than  do  the  ARS,  i.e., 

ARS  are  smoother.  Above  all,  both  differences  could  be  affected  by  the  fre¬ 
quency  increment  and  the  smoothing  technique. 

The  average  and  average  plus  one  standard  deviation  power  density  spectra 
functions  for  different  site  conditions  have  been  normalized  to  a  unit  area. 
Based  on  these  normalized  PSD  standard  spectra  and  amplification  scaling 
curves,  which  are  derived  from  the  average  power  l2  ,  the  area  under  G(uj)  , 
or  its  square  root,  the  rms  (root-mean-square)  value  (X)  developed  in  this 
study, ;a  design  earthquake  PSD  spectrum  for  any  earthquake  magnitude  and  dis¬ 
tance  can  be  generated.  These  scaling  curves,  i.e.,  the  correlation  curves 
of  peak  ground  acceleration  or  peak  velocity  versus  rms  value,  and  rms  value 
versus  Modified  Mercalli  Intensity,  can  also  serve  as  earthquake  engineering 
intensity  scales  in  a  quantitative  manner. 
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PREFACE 
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gram,  Project  4A161102AT22 ,  Work  Unit  00296,  which  was  monitored  for  the 
Office,  Chief  of  Engineers,  U.  S.  Army,  by  Mr.  A.  F.  Muller. 

This  report  was  prepared  by  Mr.  Frank  K.  Chang  of  the  Earthquake 
Engineering  and  Geophysics  Division  (EE&GD) ,  Geotechnical  Laboratory 
(GL) ,  under  the  general  direction  of  Dr.  A.  G.  Franklin,  Principal 
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McAnear,  Acting  Chief,  GL. 
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visiting  Research  Engineer  from  the  Port  and  Harbour  Research  Institute 
(PHRI) ,  Japan,  who  provided  the  digitized  Japanese  strong-motion 
accelerograms  of  1963  to  1975  recorded  by  the  PHRI  and  the  geological 
site  conditions  associated  with  each  record;  Professor  M.  Shinozuka, 
Columbia  University,  who  provided  the  computer  program  used  to  calculate 
the  one-sided  PSD  functions;  Professor  0.  W.  Nuttli,  Saint  Louis  Uni¬ 
versity,  who  made  some  comments  for  the  first  draft  of  this  report; 
Professor  E.  H.  Vanmarcke  and  Dr.  Shih-Sheng  (Paul)  Lai,  Massachusetts 
Institute  of  Technology;  Professor  H.  C.  Shah,  Stanford  University;  Dr. 

S.  C.  Liu,  National  Science  Foundation;  and  Drs.  P.  F.  Hadala,  W.  F. 
Marcuson,  and  A.  G.  Franklin,  WES,  who  made  the  critical  reviews  for 
this  final  report. 
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SITE  EFFECTS  ON  POWER  SPECTRAL  DENSITIES 
AND  SCALING  FACTORS 


PART  I:  INTRODUCTION 


Background 


1.  Three  ways  of  expressing  the  characteristics  of  earthquake 
ground  motion  in  the  frequency  domain  are  to  compute  the  peak  response 
spectrum,  Fourier  spectra,  and  power  spectral  density  (PSD).  The  re¬ 
sponse  spectrum  represents  the  peak  response  of  a  linear  1-degree-of- 
freedom  system  to  the  ground  motion  and  provides  a  convenient  method  of 
obtaining  a  preliminary  analysis  of  certain  structural  dynamics  prob¬ 
lems.  The  PSD  is  related  to  the  Fourier  amplitude  spectrum  of  the 
ground  motion.  In  fact,  the  three  are  closely  related. 

2.  Both  response  spectrum  and  PSD  functions  have  been  widely  used 
to  describe  earthquake  ground  motions  in  the  frequency  domain.  Pereira, 
Oliveira,  and  Duarte  (1977),  Vanmarcke  and  Cornell  (1972),  and  Vanmarcke 
and  Gasparini  (1977)  have  shown  that  the  PSD  function,  the  acceleration¬ 
time  plots,  and  the  response  spectra  are  all  interrelated.  Thus,  if 

a  design  level  response  spectrum  is  specified,  then  an  equivalent  PSD 
function  and  a  set  of  corresponding  acceleration-time  curves  that  are 
consistent  with  the  response  spectrum  can  be  generated.  However, 
neither  the  PSD  nor  the  response  spectrum  can  uniquely  determine  the 
acceleration-time  curve,  because  both  are  incomplete  representations  of 
the  ground  motion,  lacking  phase  information. 

3.  Seed  and  Idriss  (1971),  Tezcan  (1971),  and  many  others  have 
confirmed  that  damage  to  buildings  during  past  earthquakes  has  been 
closely  associated  with  the  vibrational  characteristics  of  the  underly¬ 
ing  soils.  Duke  and  Hradilek  (1977)  and  Berrill  (1977)  utilized  Fourier 
spectra  to  study  the  effects  of  local  site  conditions  on  ground  motions 
recorded  in  the  9  February  1971  San  Fernando  earthquake  and  its  after¬ 
shocks.  However,  no  positive  correlations  were  found.  Hudson  (1972) 


and  Hudson  and  Udwadia  (1974)  made  a  comparison  of  measured  strong 
ground  motions  in  the  form  of  plots  of  Fourier  spectra  at  selected 
sites  and  found  that  various  governing  factors  could  be  individually 
studied;  and  in  this  way,  it  was  shown  that  local  variations  are  largely 
governed  by  (a)  source  mechanisms,  (b)  propagation  path,  and  (c)  local 
geology. 

4.  Seed,  Ugas,  and  Lysmer  (1976)  presented  the  results  of  a  sta¬ 
tistical  analysis  of  the  response  spectral  shapes  of  104  ground-motion 
accelerograms  obtained  from  23  earthquakes,  mostly  in  the  western  United 
States.  The  analysis  shows  clear  differences  in  spectral  shapes  for 
different  soil  and  geological  conditions.  Chang  and  Krinitzsky  (1977) 
studied  the  duration  and  spectral  content  of  strong-motion  records  from 
the  western  United  States  according  to  site  conditions  and  found  that 
the  predominant  frequencies  are  in  the  range  of  0.1  to  6.67  Hz  and  the 
spectral  shape  depends  on  the  source  spectrum  function  (magnitude) , 
distance,  and  geological  conditions. 

5.  An  extensive  study  was  made  at  the  strong-motion  station  sites 
in  Ferndale  and  El  Centro,  California,  by  Shannon  and  Wilson,  Inc./ 
Agbabian  Associates  (1976).  Both  sites  are  located  in  a  highly  seismic 
region  and  have  many  ground-motion  records.  The  subsurface  conditions 
at  these  sites  were  defined  by  geotechnical  investigations  that  in¬ 
cluded  boring  and  sampling  of  the  subsurface  soil  materials  followed 

by  field  and  laboratory  tests  to  define  the  index  properties  and  the 
dynamic  properties  of  the  soil  needed  for  one-dimensional  wave  propaga¬ 
tion  analyses  by  numerical  methods.  The  site-response  analysis  (SHAKE 
code),  site-matched  records,  Seed-Ugas-Lysmer  site  dependent  spectrum. 
Nuclear  Commission  Regulatory  Guide  1.60  spectrum,  and  spectra  from 
the  measured  records  that  correspond  to  the  "design"  earthquake  event 
for  the  site  were  studied.  No  single  method  yielded  the  best  vibratory 
motion  criteria  for  both  sites.  The  methods  described  above  have  limita¬ 
tions  that  either  are  related  to  the  simplified  models  used  for  site- 
response  wave  propagation  analyses  or  are  a  direct  result  of  the  limita¬ 
tions  in  the  current  library  of  strong-motion  records. 

6.  Considering  the  earthquake  ground  motion  to  be  random  in 


nature,  Arnold  (1975)  and  Arnold  and  Vanmarcke  (1977)  studied  the  in¬ 
fluence  of  site  azimuth  relative  to  source  fault  orientation  and  local 
soil  conditions  on  earthquake  ground-motion  spectra  for  the  San  Fernando, 
California,  earthquake  of  9  February  1971  using  PSD  functions.  The  re¬ 
sult  showed  that  local  soil  conditions  and  site  azimuth,  as  well  as 
epicentral  distance,  can  have  a  significant  effect  on  both  the  intensity 
and  the  frequency  content  of  ground  motions.  This  study  demonstrated 
the  potential  value  of  PSD  functions  as  a  tool  for  comparing  and  study¬ 
ing  variations  in  ground-motion  characteristics  and  also  showed  the 
great  utility  of  PSD  functions  as  input  to  random  vibration  analyses  of 
structural  response. 


Purpose  and  Scope 

7.  The  main  purpose  of  this  study  is  to  analyze  the  site  depen¬ 
dence  of  PSD  shapes  in  the  frequency  range  of  0  to  10  Hz  and  to  consider 
the  application  of  PSD  shapes  in  seismic  design.  Knowledge  of  the  in¬ 
fluence  of  site  conditions  on  the  characteristics  of  earthquake  ground 
motions  and  their  spectral  shapes  is  necessary  for  earthquake-resistant 
design  and  analysis  of  structures  such  as  earth,  rockfill,  or  concrete 
dams,  large  buildings,  nuclear  power  plants,  and  other  military  or  civil 
facilities. 

8.  Factors  affecting  the  ground  motion  at  a  particular  site  in¬ 
clude  the  source  mechanism  (nature  of  fault  movement  and  magnitude  of 
energy  release),  propagation  path  characteristics,  the  direction  of  the 
site  relative  to  the  fault  rupture,  and  local  geological  and  soil  con¬ 
ditions.  This  study,  however,  deals  only  with  the  influence  of  local 
geological  and  soil  conditions  on  ground  motion. 
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PART  II:  POWER  SPECTRAL  DENSITY 


9.  In  the  application  of  the  random  vibration  theory  of  linear 
systems  for  evaluation  of  the  effects  of  variations  in  ground-motion 
characteristics  on  structural  response  to  earthquake  excitation,  the 
ground  motion  may  be  described  in  the  form  of  a  PSD  function.  The  PSD 
function  G(w)  is  defined  as  a  measure  of  the  ground-motion  power  or 
energy  per  unit  time  as  a  function  of  frequency  w  (Figure  1).  Usually, 
estimates  of  the  PSD  are  obtained  from  the  squared  amplitudes  of  the 
Fourier  transform,  or  the  squared  Fourier  amplitude  spectrum.  In  Fig¬ 
ure  1,  is  the  amplitude  of  sinusoidal  waves. 


Figure  1.  The  PSD  function  G(w) 


Fourier  Transform 


10.  Generally,  the  given  earthquake  ground-motion  function,  such 
as  an  acceleration  versus  time  plot,  can  be  represented  in  the  time  do¬ 
main  as  a(t)  and  in  the  frequency  domain  as  F(w)  : 


6 


dt  ,  i  =  v^T 


(1) 


F(oj)  =  /  a(t)e 


o 

-*/ 


F(w)e^at:  dm 


Equations  1  and  2  are  called  the  Fourier  transform  pair;  i.e.,  F(w)  >  is 

the  Fourier  integral  or  Fourier  transform  of  a(t)  ,  and  a(t)  is  the 

inverse  Fourier  transform  of  F(w)  .  The  symbols  t  and  t  denote 

the  instant  in  time  and  total  duration,  and  w  ,  w,  ,  and  w  repre- 

1  o 

sent  the  frequency,  lower  frequency  bound  =  2tt / 1  ,  and  maximum  fre¬ 

quency  (in  radians  per  second),  respectively.  For  practical  purposes, 
can  usually  be  taken  as  zero. 

Total  Intensity,  Average  Power,  and  PSD 

11.  By  Parseval's  theorem,  the  relation  between  the  energy  of 
the  motion  as  expressed  in  the  time  domain  and  in  the  frequency  domain 
can  be  represented  in  the  following  equations.  For  a  nonperiodic  func¬ 
tion,  the  total  energy  or  intensity  I  delivered  by  the  source  is 
given  by 


o  o 

Lo  -  J  |a(t>  |2  dt  -  ij 


F(oj)  dto 


The  mean-square  average  value  is  expressed 


o 

U 


I  a(t) I  dt  = 


o 


|  F (u>)  |  do)  -  /  G(oj)  da 


7 


r 


in  which  G(u>)  is  the  energy  per  unit  time  (power),  or  the  power  spec¬ 
tral  density  of  the  function  a(t)  .  The  integrand  on  the  right  side  of 
Equation  4  can  be  written  as 


G(W)  =  i  i-  j  F(ut) | 2 
o 

if  Si  (duration  of  strong  motion)  is  substituted  for  t^  ,  and 


(5) 


o 

1  r  i 

'j " 


(u>)|  div  =  S  /  G(oj)  duJ 
o 


o 

/ 


The  quantity  j  F (oj ) | “  is  called  the  energy  spectrum  or  energy  spectral 

density  function  of  a(t)  (Hsu  1967).  The  left  side  of  Equation  4 

2 

gives  the  statistical  average  power  A  of  the  function  a(t)  over 


the  total  duration  of  the  motion  t 
,  2 


o 


From  Equation  4,  also  note  that 


is  equal  to  the  area  under  the  curve  of  G(<.>)  in  Figure  1.  There- 
fore,  the  average  power  can  be  written  as 


o 

/ 


A^  =  /  G(w)  dw 


(6) 


i 


Strong-Motion  Duration  and  Normalized  PSD 


12.  From  the  relation  of  Equations  3,  4,  5,  and  6,  the  total 
intensity  I  can  be  expressed  as 


I  =  S  A 
o  o  o 


(7a) 


S  =  -a 

°  A2 

O 


(7b) 


A  peak  factor  r  may  be  defined  as 


r  = 


max 


(8a) 


A  =  — 
o  r 


(8b) 


Therefore,  by  substitution 


,  I 

2  o 


max 


(9) 


where  r  ,  which  is  a  dimensionless  parameter,  can  be  determined  empiri¬ 
cally.  This  equation  gives  the  strong-motion  duration  Sq  ,  defined  by 
Vanmarcke  and  Lai  (1977) ,  which  is  necessarily  smaller  than  the  total 


duration  t 


The  PSD  values  in  this  report  are  computed  in  terms  of 


tQ  ,  which  yields  smaller  PSD  values. 

13.  A  more  effective  way  for  dealing  with  the  frequency  content 
of  ground  motion  is  through  the  normalized  spectral  density  function 
G*(w) : 


G*(u)  =  -4  G(u)  (10) 

A 

o 

If  G*(to)  is  an  individual  normalized  PSD  function,  the  statistical 
mean  PSD  curve  will  be 


n 

G* (to)  =  -  y  G*(w)  ,  i  =  1,2,... n  (11) 

i=l 

In  practice,  curves  of  G*(oo)  computed  from  suites  of  actual  earthquake 
records  show  large  and  irregular  fluctuations.  To  isolate  the  system¬ 
atic  component  from  the  random  variations,  frequency  smoothing  with  a 
"Hanning"  window  (Blackman  and  Tukey  1958)  should  be  used.  A  detailed 


procedure  for  the  frequency  smoothing  will  be  presented  in  Part  III. 


Average  acceleration: 


X 

o 


t 

o 

1/2 

0) 

o 

r 

= 

~  f  |a(t) |2  dt 

X  = 

n 

f  G(o j)  dw 

Jo 

1/2 


(15a,  b) 


Strong-motion  duration: 

I 


o 


(16) 


15.  This  study  includes  investigation  of  (a)  the  site  dependence 
of  the  PSD's  or  the  power  spectral  shapes  in  the  frequency  range  of 
0  to  10  Hz  for  four  general  types  of  site  conditions,  which  are  classi¬ 
fied  as  rock,  stiff  soil,  deep  cohesionless  soil,  and  soft  soil;  (b)  the 
statistical  characteristics  of  the  ground  motions;  (c)  the  relations  be¬ 
tween  the  power  spectrum,  average  acceleration,  average  power  (scaling 
factor),  duration,  and  total  intensity  (Arias);  and  (d)  the  development 
of  average  acceleration  as  an  alternative  earthquake  engineering  inten¬ 
sity  scale.  The  average  acceleration  not  only  describes  the  intensity 

of  ground  motion  for  input  to  structural  design  but  also  serves  as  the 
2 

scaling  factor  (Xq)  of  the  normalized  standard  PSD  spectra.  In  this  re¬ 
port,  duration  will  be  considered  in  a  general  sense;  the  strong-motion 
duration  will  not  be  defined. 
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PART  III:  DATA  SELECTION  AND  DATA  PROCESSING 


Data  Selection 

16.  A  total  of  421  horizontal  ground  accelerograms  from  89  earth¬ 

quakes,  mostly  in  the  western  United  States  and  Japan  (will)  a  few  rec¬ 
ords  from  Russia,  Rumania,  and  India),  were  selected  for  this  analysis. 
Based  on  the  site  classifications  of  Seed  and  Idriss  (1971)  and  Seed, 
Ugas ,  and  Lysmer  (1976),  these  records  have  been  divided  into  four 
groups:  (a)  56  records  for  rock  sites,  (b)  131  records  for  stiff  soil 

sites  (depth  <150  ft),  (c)  120  records  for  deep  cohesionless  soil  sites 
(depth  >250  ft),  and  (d)  114  records  for  soft  to  medium  clays  with  asso¬ 
ciated  strata  of  sands  or  gravels.  One  hundred  seventy-three  of  the 
421  records  were  obtained  from  California  Institute  of  Technology  (CIT) 
Volume  II-Corrected  Accelerograms  (1971-75),  and  220  uncorrected  records 
were  provided  by  the  Port  and  Harbour  Research  Institute  (PHRI),  Japan. 
The  digitized  Gazli  (USSR)  and  Bucharest  (Rumania)  records  were  provided 
by  Dr.  A.  G.  Brady,  U.  S.  Geological  Survey.  All  421  corrected  and  un¬ 
corrected  records  were  adjusted  to  zero  mean  before  processing  the  PSD. 

Definition  of  Average  Power  and  Average  Acceleration 

17.  The  main  approach  used  in  this  study  was  to  determine  the 
normalized  mean  and  the  mean  plus  one  standard  deviation  PSD  shape 
(NPSD)  for  each  group,  and  the  average  acceleration  1  using  Equa¬ 
tion  15b,  for  each  raw  record.  Figure  1  shows  G(o i)  ,  whose  value 

2  2 

at  w.  is  equal  to  A./2Aui  ,  so  that  A  is  actually  the  power  or 

1  1  O  " 

energy  density  in  an  accelerogram  for  a  finite  frequency  band  (0  £  f 
£  10  Hz  in  this  study),  and  A ^  is  the  amplitude  of  the  ith  wave  compo¬ 
nent  in  centimetres  per  second  squared.  The  total  or  average  power  will 
become  equal  to  the  area  under  the  continuous  curve  C(ui)  . 

Spectra,  Smoot hing 

18.  The  statistical  (mean  and  mean  plus  one  standard  deviation) 
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m 


NPSD  function  shapes  for  each  group  appear  very  irregular.  Therefore, 
a  spectral  or  frequency  smoothing  technique  has  to  be  employed  to  elimi¬ 
nate  random,  or  nonsystematic ,  fluctuations  of  the  NPSD  curve.  A  final 
smooth  estimate  of  the  NPSD  may  now  be  formed  by  further  frequency 
smoothing  with  a  procedure  called  the  "Hanning  window"  ly  Blackman  and 

Tukey  (1958)  and  Bendat  and  Piersol  (1971).  Let  G  and  G  denote  a 

K.  K 

raw  and  smooth  estimate  at  harmonic  K  ,  where  K  =  0,l,2,...m  ;  then 

G  =  0.5G  +  0.5G, 
o  o  1 

Gr  =  0.25Gk_1  +  0.5Gr  +  0.25Gr+1  K  =  l,2,...m  +  1  (17) 

G  =  0.5G  .  +  0.5G 
m  m~l  m 

Statistical  Errors 


19.  The  descriptive  properties  of  a  random  variable  cannot  be 
precisely  determined  from  sample  data.  Only  estimates  of  the  parameters 
of  interest  can  be  obtained  from  a  finite  sample  of  observations.  The 
accuracy  of  parameter  estimates  based  upon  sample  values  can  be  described 
by  a  mean  square  error  defined  as 


e[(<!>  -  <5>)2]  =  E  |  [  <t>  -  E($)3  2| 


(18) 


where  <t>  is  an  estimator  for  i>  .  The  first  term  on  the  right  side  of 

Equation  18  is  the  variance  Var  (<J>)  ,  which  describes  the  random  por- 

2  " 

tion  of  the  error;  the  second  term  is  the  square  of  a  bias  b  (<f)  , 
which  describes  the  systematic  portion  of  error.  Therefore,  the  mean 
square  error  is  the  sum  of  two  terms: 


=  Var  (4>)  +  b2  (4>) 


and  the  rms  error  is 


o2  (<J>)  +  b2  (40 


(19) 


(20) 
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where  Var  (4>)  =  a  (4>)  and  a(4>)  equals  the  standard  error  or  random 

error.  Bendat  and  Piersol  (1971)  give  the  simple  relationship  between 

the  random  error  E  and  the  smoothing  times  N,  as  E  =  I/^nT  • 
r  d  r  d 

Thus,  for  increasing  numbers  of  smoothing  times,  the  corresponding  values 

of  E  are  as  follows: 
r 


0.045 


10,000 


The  four  normalized  mean  and  mean  plus  one  standard  deviation  PSD  curves 
in  this  study  have  been  smoothed  500  times,  so  the  random  error 
should  be  less  than  5  percent.  Figure  2  shows  the  effect  of  the  smooth¬ 
ing  process  on  a  typical  PSD. 


Record  Length,  Increment  Frequency,  and  NPSD  Function 

20.  The  record  length  and  spectral  content  (spectral  amplitude 

and  frequency  range)  are  two  basic  elements  for  controlling  the  spectral 

intensity.  The  incremental  frequency  Af  ,  used  in  the  PSD  computer 

program,  depends  on  the  total  record  length.  Since  it  is  necessary  to 

use  the  same  value  for  Af  throughout,  all  accelerograms  have  been 

13 

processed  to  give  them  a  duration  of  163.82  sec  or  8192  (2  )  digital 

points  with  an  equal  time  interval  At  of  0.02  sec,  which  gives 
Af  =  0.006104  Hz  .  Outside  the  time  of  the  actual  record,  the  ampli¬ 
tudes  at  extended  times  were  set  to  zero.  In  this  study,  the  PSD  func¬ 
tion  G(f)  has  been  defined  to  include  only  the  frequency  range  of 
0  to  10  Hz  so  that 


10  Hz 


G (f )  df 


Since  Af  is  0.006104  Hz,  there  are  1640  points  in  the  raw  PSD  function 
for  each  accelerogram.  The  NPSD  function  is  defined  as  the  PSD  amplitude 
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divided  by  the  area  under  the  power  density  versus  frequency  curve. 

This  area  may  also  he  called  the  spectral  intensity. 

Spectral  Frequency  Range 

21.  Figures  2a,  b,  and  c  show  the  PSD  spectra  of  N-S,  E-W,  and 
vertical  components,  respectively,  of  the  El  Centro  earthquake  of 

18  May  1940.  The  frequency  range  shown  is  between  0  and  24  Hz.  While 
the  energy  of  the  PSD  in  the  vertical  component  spreads  between  0  and 
22  Hz,  the  energy  of  the  PSD  for  the  two  horizontal  components  is  con¬ 
centrated  in  the  range  of  0  to  10  Hz .  For  this  record,  the  PSD  fre¬ 
quency  range  of  0  to  10  Hz  is  adequate  for  the  description  of  the 
spectra  of  the  horizontal  components. 

Computer  Procedures  for  Generating  PSD 

22.  The  computer  program  used  to  calculate  one-sided  PSD  function 
was.  provided  by  Professor  M.  Shinozuka,  Columbia  University,  and  modi¬ 
fied  for  the  Honeywell  635  Computer  by  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  Automatic  Data  Processing  (WES  ADP)  Center.  The 
procedures  for  generating  PSD  are  as  follows: 

a.  Read  in  accelerogram. 

In  Scale  accelerogram*  so  that  accelerations  are  in  centi¬ 
metres  per  second  squared. 

£.  Adjust  accelerogram*  by  interpolation  so  that  the  At 
between  time  points  is  0.02  sec. 

13 

cL  Extend  accelerogram  to  (2  -  1)  x  0.02  =  163.82  sec  by 

adding  trailing  zeros  to  the  acceleration  record  (so  far, 
no  accelerogram  is  longer  than  163.82  sec). 

e_.  Calculate  the  statistics  (mean,  standard  deviation,  etc.) 
for  the  extended  accelerogram  and  adjust  it  to  a  zero 
mean. 

f_.  Calculate  the  PSD  of  the  extended,  zero-mean  accelerogram 
(up  to  a  frequency  of  10  Hz). 


*  If  necessary. 
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Calculate  the  area  under  the  PSD  curve  and  normalize  the 
curve  (NPSD)  (i.e.,  divide  each  point  of  the  PSD  by  the 
area  under  the  PSD  curve). 

Smooth  the  PSD  curve  or  NPSD  curve  using  the  "Hanning" 
process. 


PART  IV:  DATA  PRESENTATION 


23.  A  total  of  421  horizontal  accelerograms  were  used  to  estimate 
the  raw  power  spectral  density  (RPSD)  function  and  its  variance  or  RPSD 
intensity.  These  calculated  PSD  functions,  normalized  to  unit  area, 
will  be  used  to  represent  the  frequency  content  and  characteristics  of 
ground  motions.  The  final  mean  normalized  PSD  (MNPSD)  and  MNPSD  plus 
one  standard  deviation  curves  for  four  site  groups  were  plotted. 

24.  Depending  on  their  recording  site  conditions,  the  421  horizon¬ 
tal  records  were  grouped  in  Tables  1  to  4: 

Table  1-56  records  for  rock  sites 
Table  2  -  131  records  for  stiff  soil  sites 

Table  3  -  120  records  for  deep  cohesionless  soil  sites 

Table  4  -  114  records  for  soft  to  medium  clay  and  sand 

The  headings  for  Columns  1-9  in  each  table  are  self-explanatory.  The 

duration  in  Column  10  is  arbitrarily  estimated.  Column  11  gives  the 

2  2 

base  PSD  average  power  A  (note  that  A  =  1^/164),  which  is  the  area 
under  the  RPSD  function  curve  for  an  extended  record  length  of  163.82 
sec,  or  164  sec  for  simplicity.  Column  12  gives  A  ,  the  average  accel¬ 
eration,  or  the  square  root  of  the  corresponding  value  in  Column  11. 

2  9 

Column  13  is  the  factor  for  conversion  of  A  (Column  11)  to  A"  ,  the 

o 

area  under  the  PSD  function  for  the  actual  record  length  (Column  14). 

This  conversion  factor  in  Column  13  is  the  ratio  of  the  extended  time  of 

163.82  sec  to  the  time  of  the  total  record  length  t  or  the  arbi- 

o 

trarily  selected  duration  (Column  10);  i.e. ,  Column  11  *  Column  13 

2 

=  Column  14,  which  is  the  raw  average  power  A  .  Column  15  gives  X  , 

o  o 

the  average  acceleration  for  the  corresponding  Column  14.  Column  16 

equals  Column  14  multiplied  by  0.875  and  is  the  final  corrected  area,  or 
2 

average  power  Ag  under  the  RPSD  function  curve.  The  constant  0.875  is 
the  amount  by  which  the  power  spectrum  estimates  should  be  multiplied  so 
as  to  obtain  the  correct  scale  factor.  An  explanation  for  this  correc¬ 
tion  is  presented  by  Bendat  and  Piersol  (1971,  p.  323).  Column  17  gives 
the  average  acceleration  Ag  or  the  square  root  of  the  area  under  the 
estimated  PSD  function  curve.  Actually,  the  difference  between  Aq  and 


A  is  very  small,  so  it  can  be  considered  that  A  s  A  .  For  practical 
s  ^  o  s  r 

purposes,  the  raw  average  power  A^  may  be  accepted  as  the  scaling  fac¬ 
tor  for  the  standard  (mean  or  mean  plus  one  standard  deviation)  NPSD 
spectrum. 

25.  Most  of  the  records  used  in  this  study  were  obtained  at  sites 
in  the  western  part  of  the  United  States  or  Japan.  A  few  other  signifi¬ 
cant  strong-motion  records,  such  as  those  of  the  Koyna,  India  (1967), 
Gazli,  USSR  (1976),  and  Bucharest,  Rumania  (1977)  earthquakes,  were  also 
included. 

26.  The  RPSD  average  power  values  of  the  horizontal  components  for 
the  earthquake  accelerograms  of  El  Centro  (1934  and  1940),  Taft  (1952), 
an'1  Olympia  (1949)  have  been  calculated  for  durations  of  25  or  30  sec  as 
indicated  in  Table  5.  The  average  power  of  uncorrected  versions  of  these 
records  (for  the  same  durations)  was  calculated  by  Ravara  (1965)  and 
should  be  different  from  that  calculated  in  this  study.  For  the  purpose 
of  comparison  and  evaluation  of  accuracy,  Ravara’ s  values  are  also  listed 
in  Row  1.  The  numerical  values  in  Row  2  are  directly  calculated  from 
the  CIT  corrected  records.  Average  power  values  adjusted  from  these 
calculations  for  a  duration  of  163.82  sec  to  the  lengths  of  the  particu¬ 
lar  records  are  shown  in  Row  3.  In  comparing  the  average  power  in  Row  2 

and  Row  4  of  Table  5,  there  is  a  close  agreement.  The  average  power 

2 

values  in  Row  3  estimated  from  the  base  average  power  A  are  higher 
than  those  in  Row  2,  which  were  directly  calculated  from  the  actual 
(shorter)  duration  records.  Thus,  the  amount  of  increased  average  power 
(~  a  constant  factor)  caused  by  adding  zeros  indicates  that  the  final 
correction  is  needed. 

27.  Table  5  also  shows  that  the  average  power  estimated  from  the 
baseline  uncorrected  accelerograms  by  Ravara  (1965)  could  have  18  to  66 
percent  error  in  comparison  with  CIT  corrected  data.  The  error  for  the 
extended  163.82-sec  record  duration  in  this  study  is  in  the  range  of 

8  to  16  percent.  The  average  error  is  about  12  percent,  which  is  in 
agreement  with  the  correction  factor  of  0.875  (Bendat  and  Piersol  1971). 
After  the  correction  factor  of  0.875  is  applied,  the  error  is  reduced 
from  0.3  to  5.5  percent. 


2] 


28.  In  conclusion,  the  base  average  power  X  of  421  records 
grouped  in  Column  11  of  Tables  1-4,  which  are  presented  in  a  convenient 
way,  can  be  employed  to  estimate  the  average  power  for  any  strong-motion 
duration  of  an  individual  accelerogram  in  the  lists  of  all  four  tables  as 
long  as  the  selected  duration  is  less  than  163.82  seconds.  Fortunately, 
not  one  of  the  421  records  exceeds  the  duration  of  163.82  seconds. 

29.  Table  6  is  a  comparison  of  average  power  calculated  by  the 

method  of  Vanmarcke  and  Lai  (1977)  and  Vanmarcke  (1979)  and  the  method 

of  this  study  with  the  same  strong-motion  durations  corresponding  to  the 

2 

same  records.  By  comparison  of  the  average  power  Xq  and  the  average 

acceleration  X  for  the  same  record  in  Table  6,  it  follows  that  the 
o 

values  calculated  by  Vanmarcke  and  Lai  (1977)  and  in  this  study  are  in 

2 

excellent  agreement,  even  though  Vanmarcke  and  Lai's  values  of  X^  were 
directly  derived  from  the  time  domain  and  the  values  in  this  study  are 
calculated  in  the  frequency  domain.  This  comparison  is  verifying  not 
only  the  processing  techniques  but  also  the  theoretical  background. 
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PART  V:  DATA  ANALYSES 


Analysis  of  Site-Dependent  PSD  Spectral  Shape 

30.  The  NPSD  functions  for  the  horizontal  accelerograms  in  each 
group  (Table  1-4)  were  first  determined  and  then  were  analyzed  statis¬ 
tically  to  obtain  the  average  NPSD  spectra  and  the  average  plus  one 
standard  deviation  NPSD  spectra  (about  84  percentile).  Figures  3-6 
present  these  mean  and  mean  plus  one  standard  deviation  NPSD  spectra 
for  the  four  different  site  conditions.  The  mean  NPSD  spectra  for  the 
different  site  conditions  are  compared  in  Figure  7,  and  the  mean  plus 
one  standard  deviation  NPSD  spectra  are  compared  in  Figure  8.  Both 
NPSD  spectra  in  Figures  7  and  8  are  smoothed  500  times. 

31.  It  is  clear  that  the  differences  in  PSD  spectral  shapes  de¬ 

pend  on  the  site  conditions.  In  particular,  two  categories  can  be 
distinguished:  sites  having  soft  to  medium  clays  and  sands  or  deep 

cohesionless  soils  (>250  ft)  are  similar  and  form  one  category,  the 
soft  group;  stiff  soil  and  rock  sites  form  another  category,  the  hard 
group.  A  dividing  line  on  the  frequency  axis  appears  at  2.5  Hz  (0.4-sec 
period).  In  the  frequency  range  below  2.5  Hz,  spectral  amplifications 
for  the  soft  group  are  much  higher  than  those  for  the  hard  group;  in 
the  frequency  range  above  2.5  Hz,  spectral  amplifications  for  the  hard 
group  are  higher  than  those  for  the  soft  group.  In  the  soft  group, 

the  energy  peaks  for  the  soft  to  medium  clays  and  sands  and  the  deep 
cohesionless  soils  both  occur  at  about  the  same  frequency  of  1  Hz,  but 
the  amplifications  are  slightly  different  (0.4  and  0.35,  respectively). 

32.  The  average  NPSD  spectrum  of  the  deep  cohesionless  soil  sites 
has  a  large  hump  at  2.8  Hz  (0.36-sec  period),  but  the  spectrum  of  the 
soft  to  medium  clay  and  sand  sites  does  not.  The  large  amplitude  at 

0  Hz  is  believed  to  be  caused  by  a  digitization  error,  particularly  that 
due  to  the  uncorrected  Japanese  strong-motion  data  in  the  soil  site 
group.  For  the  hard  site  group,  in  the  frequency  range  below  2.5  Hz, 
the  spectral  amplitude  for  the  stiff  soil  is  higher  than  that  for  the 
rock;  but  in  the  frequency  range  above  2.5  Hz,  the  spectral  amplitude 
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Figure  3.  Mean  and  mean  plus  one  standard  deviation  (o) 
NPSD  curves  of  rock  sites,  raw  (left)  and  smoothed  500 

times  (right) 
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Figure  4.  Mean  and  mean  plus  one  standard  deviation  (o) 
NPSD  curves  of  stiff  soil  sites,  raw  (left)  and  smoothed 

500  times  (right) 
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Figure  7.  Comparison  of  mean  NPSD  curves  of  four  soil  groups 
raw  (left)  and  smoothed  500  times  (right) 
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Figure  8.  Comparison  of  mean  plus  one  standard  deviation 
NPSD  curves  of  four  soil  groups,  raw  (left)  and  smoothed 

500  times  (right) 
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for  the  rock  is  slightly  higher  than  that  for  the  stiff  soil.  The 
largest  energy  peaks  for  the  rock  sites  and  the  stiff  soil  sites  are 
at  2.75  Hz  (0.36  sec)  and  0.8  Hz  (1.25  sec),  respectively. 

Statistical  Characteristics  of  the 
Earthquake  Ground  Motions 

33.  Based  on  the  maximum  ground  accelerations  and  average  ac¬ 
celerations  (PSD  intensities)  in  Tables  1-4,  and  the  average  NPSD  func¬ 
tion  estimates  for  the  four  site  conditions  (Figures  7  and  8),  Table  7 
lists  the  statistical  characteristics  of  the  earthquake  ground  motions. 

34.  It  is  apparent  from  Table  7  that  statistical  characteristics 
of  the  ground  motions  are  strongly  site  dependent.  The  rock  sites  pro¬ 
duce  an  average  maximum  ground  acceleration  of  about  0.20  g  for  the 
entire  suite  of  56  records.  The  rock  site  group  shows  the  highest 
maximum  acceleration  and  average  acceleration  of  the  four  site  groups. 
The  PSD  function  estimates  are  almost  uniform  over  the  peak  frequencies 
of  1.06,  2.75,  3.80,  and  5.17  Hz. 

35.  The  average  maximum  ground  accelerations  and  PSD  spectral 
intensities  (or  average  accelerations)  for  the  other  three  site  groups- 
stiff  soils,  deep  cohesionless  soils,  and  soft  to  medium  clays  and 
sands — are  relatively  close  together.  However,  the  spreads  of  the 
standard  deviation  of  maximum  accelerations  for  the  stiff  soil  and 
cohesionless  soil  groups  are  wider  than  for  the  soft  to  medium  clays 
and  sands  group.  This  large  spread  is  believed  to  be  caused  by  the 
different  magnitudes  of  earthquakes  and  the  different  epicentral  dis¬ 
tances.  The  group  of  accelerograms  for  soft  to  medium  clays  and  sands 
shows  relatively  low  ground  acceleration  but  the  highest  PSD  function 
estimates  at  the  frequency  of  1  Hz  among  the  four  groups.  One  second 
(1  Hz)  is  probably  near  the  predominant  natural  period  of  sites  on  soft 
to  medium  clays  and  sands.  It  seems  that  the  acceleration  and  the  PSD 
spectral  intensity  at  1  Hz  are  roughly  in  proportion  and  inverse  proper 
tion,  respectively,  to  the  degree  of  stiffness  of  the  site  material. 

In  conclusion,  the  average  acceleration  or  the  mean  spectral  intensity 
is  site  dependent. 
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Maximum  Ground  Acceleration  and  Average  Acceleration 

36.  Figures  9-16  show  plots  of  maximum  ground  accelerations  a 

max 

against  base  average  accelerations  A  (i.e.,  the  rms  value  of  the  aver¬ 
age  power  for  extended  durations  of  163.82  sec)  and  against  the  average 

acceleration  A  (for  selected  or  actual  record  durations)  for  each  of 
s 

the  four  site  condition  groups.  These  figures  indicate  a  common  linear 
trend  for  all  four  groups.  This  approximately  linear  relationship  may 
provide  a  basis  for  predicting  strong  earthquake  ground  motions  for  engi¬ 
neering  design.  The  data  points  show  greater  scatter  on  the  plots  of 

a  versus  A  than  on  the  plots  of  a  versus  A  .  However,  mean 
max  max  s 

lines  for  both  kinds  of  plot  are  parallel,  probably  because  of  the  close 

relationship  between  A  and  Ag  .  Also,  the  data  points  for  soft  sites 

show  a  wider  spread  than  those  for  hard  sites. 

36.  It  is  worthwhile  to  note  that  there  is  a  strong  correlation 

(Figures  17  and  18)  between  a  and  A  ,  which  were  derived  from  the 
2  max  ° 

values  of  I  =  S  A  of  140  strong-motion  records  in  Vanmarcke  and  Lai 
o  o  o 

(1977,  1980).  The  quantity  I  is  the  total  motion  energy  at  constant 

2  ° 

power  A^  over  the  strong-motion  duration  Sq  .  Twenty-two  of  the  140 

records  are  for  rock  sites,  and  the  rest  are  for  soil  sites.  At  the 

same  time,  the  two  mean  lines  for  rock  sites  and  soil  sites  are  almost 

identical,  thus  indicating  that  the  linear  relationship  between  a 

°  max 

and  A  is  independent  of  the  site  conditions.  The  mean  lines  of  a 

o  max 

versus  A^  calculated  from  the  strong-motion  data  of  Vanmarcke  and  Lai 
(1977)  are  also  plotted  in  Figures  10,  12,  14,  and  16  and  lie  to  the 
right  of  the  data  for  this  study.  Since  A^  is  inversely  proportional 
to  Sq  ,  Ao  is  minimum  when  the  whole  record  length  is  considered. 

The  durations  selected  in  this  study  were  close  to  the  whole  record 
lengths;  thus,  the  calculated  average  accelerations  are  lower  than  the 
rms  accelerations  of  Vanmarke  and  Lai  (1977). 

Peak  Factor 


38.  The  peak  factor  r  is  defined  as  the  ratio  of  the  peak 
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Figure  9.  Correlation  of  maximum  ground  acceleration  (g)  versus 
"base"  average  acceleration  calculated  based  on  extended  163.82- 
sec  duration  for  rock  sites 


MAX.  GROUND  ACCELERATION,  g 


AVERAGE  ACCELERATION  OR  rms  ACCELERATION,  cm/sec2ICORRECTED' 


Figure  10.  Correlation  of  maximum  ground  acceleration  (g)  versus 
average  acceleration  and  rms  acceleration  based  on  actual  dura¬ 
tions  and  strong-motion  durations  for  rock  sites 
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Figure  12.  Correlation  of  maximum  ground  acceleration  (g)  versus 
average  accelerations  and  rms  accelerations  based  on  actual  dura¬ 
tions  and  strong-motion  durations  for  stiff  soil  sites 
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Figure  13.  Correlation  of  maximum  ground  acceleration  (g)  versus 
"base"  average  acceleration  based  on  extended  163.82-sec  duration 
for  deep  cohesionless  soil  sites 
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Figure  14.  Correlation  of  maximum  ground  acceleration  (g)  versus  average 
acceleration  and  rms  acceleration  based  on  selected  actual  durations  and 
strong-motion  durations  for  deep  cohesionless  soil  sites 
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Figure  15.  Correlation  of  maximum  ground  acceleration  (g)  versus  "base" 
average  acceleration  based  on  extended  163.82-sec  duration  for  soft  to 

medium  clay  and  sand  sites 
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MAX.  GROUND  ACCELERATION, 
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Figure  17.  Maximum  ground  acceleration  (g)  versus  rms  or  average 
acceleration  for  strong-motion  duration  case  -  rock  sites 
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Figure  18.  Maximum  ground  accelerations  (g)  versus  rms  or  average 
acceleration  for  strong-motion  duration  case  -  soil  sites 


ground  acceleration  a  to  the  average  acceleration  1  ,  or 

max  o 


The  physical  meaning  of  r  is  the  slope  of  the  lines  plotted  in  Fig¬ 
ures  9-18.  Equation  22  may  be  written  as 


X  =  — 
o  r 


If  Equation  23  is  substituted  into  Equation  16,  the  result  may  be 
expressed  as 

S  =  r2  (24) 

o  2 

a 

max 

where  r  is  a  constant  that  may  be  determined  from  Figures  9  through  16. 
The  average  peak  factors  for  rock  sites,  stiff  soil  sites,  deep  cohesion¬ 
less  soil  sites,  and  soft  soil  sites  found  in  this  study  were  5.911, 
5.422,  6.996,  and  5.695,  respectively.  Evidently ,  the  peak  factors  are 
nearly  independent  of  site  conditions.  Therefore,  an  average  peak 
factor  of  6.0  is  an  adequate  estimate  for  use  with  long  record  lengths 
such  as  those  used  in  this  report.  However,  r  is  dependent  on  the 
choice  of  record  duration.  Vanmarcke  and  Lai  (1980)  found  the  average 
peak  factor  for  140  horizontal  earthquake  records  to  be  about  2.75 
because  they  used  strong-motion  durations.  Their  simplified  definition 
of  strong-motion  duration  is 


?  T  1 

S  *  (2.75V  =7.5 

n  A  A 


However,  Equation  25  is  not  employed  in  this  studv. 

Comparison  of  Average  Acceleration  and  Leak 
Ground  Acceleration  Versus  Distance 


39.  The  average  acceleration  '  ^  in  Tables  1-3  and  the  peak 
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ground  acceleration  (PGA)  or  amax  of  the  San  Fernando  earthquake  of 
9  February  1971  are  compared  qualitatively  in  Figures  19  and  20.  The 
spread  of  the  average  acceleration  Is  as  wide  as  that  of  the  PGA,  but 
the  attenuation  of  the  average  acceleration  is  slower  than  that  of  the 
PGA.  The  average  acceleration  of  an  accelerogram  is  inversely  propor¬ 
tional  to  the  duration,  which  in  this  study  was  arbitrarily  selected. 

In  earthquake  engineering  design,  both  strong-motion  duration  and  wave 
amplitude  should  be  considered.  The  average  power  of  the  PSD  function 
and  the  average  acceleration,  the  square  root  of  the  average  power, 
possess  information  on  both  duration  and  amplitude.  The  previous  sec¬ 
tions  have  shown  good  correlation  between  the  PGA  or  a  and  the 

max 

average  acceleration  A^  .  The  average  acceleration  can  provide  an 
alternative  earthquake  engineering  intensity  scale,  describe  the  inten¬ 
sity  of  ground  motion  for  input  to  structural  design,  and  also  serve  to 

2 

compute  the  scaling  factor  A  of  the  normalized  standard  PSD  spectra. 

Potential  Uses  of  Site-Dependent  Standard 
NPSD  Spectral  Curves 

40.  In  the  previous  sections  of  Part  V,  the  four  standard  site- 

dependent  PSD  spectral  curves  at  the  common  record  length  of  163.82  sec 

have  been  established,  and  they  were  also  normalized  to  a  unit  area. 

The  relationships  between  the  shape  of  the  spectral  density  function 

and  the  duration  of  strong  ground  motion  for  individual  records  will  be 

explained  in  an  example.  It  is  very  easy  to  select  an  actual  record 

in  Tables  1-4  and  to  modify  its  duration  only  for  consistency  with  a 

specified  design  earthquake.  Use  the  N-S  component  of  the  El  Centro 

record,  18  May  1940,  as  an  illustration:  Record  No.  3  of  Table  2  shows 

2 

that  the  base  PSD  intensity  A  (average  power)  for  a  duration  of 

2  U 

163.82  sec  is  663.202  cm  /sec  .  Next,  let  the  specified  duration  for 
seismic  design  be  40  sec,  then  the  new  PSD  intensity  will  be 

,  ^  —  AA'i  909  163 . 82  991  ii  n  2.  4 

AQ  -  ooj.202  x  — —  =  2716.1  cm  /sec 

2  4 

This  value,  2716.1  cm  /sec  ,  is  also  the  scaling  factor. 
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Figure  20.  Peak  ground  acceleration  versus  distance  of  San 
Fernando  earthquake,  9  February  1971 
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41.  An  alternative  way  is  to  empty  the  approximate  linear  rela¬ 
tionships  between  a  and  X  and  between  a  and  1  (for  the 

max  max  o 

rock  site  condition.  Figures  9  and  10).  If  '.30-g  maximum  ground  accel¬ 
eration  and  15-sec  duration  are  chosen  as  the  design  ground  motion,  what 
will  the  scaling  factor  be  for  the  normal  ized  mean  PSD  curve  for  rock 

sites?  In  Figure  9,  the  base  average  acceleration  A  for  0.30  g  is 
a 

15.5  cra/sec',  and  the  conversion  factor  is  153.82/15  =  10.923;  then 

(15.5)-  x  10.923  =  2624.25  cm'~/sec^  (A  =  51.22  tm/sec") .  Thus, 

2,4  ° 

2624.25  cm  /sec  is  the  scaling  factor. 

Peak  VeJ  oc i t v  Vor sus  A verage  Acce 1  or at  ion 

42.  Figure  21  shows  the  relationship  between  the  peak  velocity 
and  the  average  acceleration.  All  140  average  acceleration  values  in 
this  figure  were  calculated  from  the  total  ground  motion  intensity  I 

o 

listed  in  the  tables  of  Vanmarcke  and  1  .  i  i  (  1977)  or  Vanmarcke  (1980), 
except  those  of  Gazli  and  Paco iron.  All  peak  velocities  are  given  hv 
Chang  (1978).  These  data  points  (4.5  M  •  6.8)  spread  over  a  wide 
band;  the  upper  line  shown  corresponds  to  magnitude  6.5  and  the  lower 
to  magnitude  5.5.  In  the  relationship  between  the  peak  acceleration 
and  the  average  acceleration  in  Figures  17  and  18,  this  particular 
feature  is  not  shown,  because  the  velocity  is  related  to  the  intensity 
or  energy  Level,  and  thus  magnitude.  The  largest  earthquake  representeu 
is  the  Kern  County  earthquake  of  1952,  for  which  the  surface-wave  magni¬ 
tude  M  was  estimated  as  7.7.  Professors  Bolt  (1978)  and  Nut t 1 i  et  al . 
s 

(1979)  found  the  local  magnitude  M  and  the  body-wave  magnitude 
to  be  7.2  and  6.8,  respectively.  Thus,  6.8  has  replaced  7.7  in 
Figure  21. 

43.  Vanmarcke  and  l.ui’s  (1077,  1980)  total  intensity  data  were 
used  for  calculating  average  acceler at  ion  because  they  have  a  unique 

definition  of  strong-motion  duration.  these  data  were  obtained  f rom 

2  ■> 

I  =  SA  where  S  is  the  strong-mot  ion  duration  and  is  the 

o  o  o  o  o 

square  of  average  acceleration,  or  average  power. 

4  1 
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Figure  21.  Correlation  of  peak  velocity  versus  average  acceleration 
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Cor relat ion  of  Average  Acceleration  and  Modified 
MercaLli  Intensity 


44.  It  will  be  of  much  benefit  to  the  engineering  community  if 

a  quantitative  earthquake  intensity  scale,  such  as  average  acceleration 
intensity,  can  be  correlated  with  the  Modified  Mercalli  Intensity  (MMI) 
(Figures  22  and  23).  '['able  8  shows  tiie  upper  bound  of  site-dependent 
rms  intensity  (square  root  of  the  sum  of  two  horizontal  average  powers) 
versus  the  MMI.  fo  obtain  these  bounds  in  the  table,  the  rock  and  stiff 
soil  sites  in  Figure  22  were  combined  as  hard  sites  and  the  deep  cohe¬ 
sionless  soil  and  soft  soil  sites  in  Figure  23  as  soft  sites.  The 
Pacoima  Dam,  Karakvr  Point,  Kovna  Dam,  and  hake  Hughes  Array  No.  12 
sites  were  located  in  the  epicentral  regions  and  near  the  faults 
(3km  £  R  *.  20km).  Certainly,  they  possessed  the  maximum  average  accele¬ 
ration  and  might  be  called  epicentral  average  accelerations.  It  seems 
from  these  limited  data  that  the  maximum  average  acceleration  at  the 

•3 

epicentral  region  might  not  he  over  550.0  cin/seo".  However,  the  power 

2 

or  the  average  acceleration  '  is  inversely  proportional  to  the 
duration,  i.e.,  the  smaller  the  duration,  the  higher  the  average 
acceleration. 

Correlation  ol  i  and  MMI 

_ _  _  o _ 

45.  Tiie  correlations  of  total  intensity  I  with  the  MMI  based 

o 

on  the  data  of  hard  sites  (Tables  1  and  2)  and  soft  sites  (Tables  3 
and  4),  are  plotted  in  Figures  24  and  25,  respectively.  The  upper 
bound  line  of  Figure  24  is  established  bv  five  earthquakes  (San  Fernando, 
C.azli,  Parktield,  Kovna,  and  Tokachi  Oki).  There  are  four  sites  lo¬ 
cated  in  the  epicentral  region  (sec  paragraph  44),  so  the  values  are 
named  as  epicentral  intensities.  The  extrapolation  from  these  values, 
the  probable  epicentral  seismic  intensities  versus  the  MMI,  may  be 
listed  as  follows: 
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Figure  24.  Probable  seismic  intensity  1  at 
epicentrnl  region  -  hard  sites  ° 
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46.  The  upper-bound  line  of  Figure  25  indicates  that  the  seismic 
intensity  at  sott  sites  is  lower  than  at  hard  sites.  Also,  the  rate  of 
attenuation  is  lower  for  soft  sites  than  for  hard  sites.  Of  course,  the 
upper-bound  seismic  intensities  for  the  hard  sites  (Figure  24)  and  the 
soft  sites  (Figure  25)  are  in  the  near  field.  The  data  under  the  upper- 
bound  line  (both  Figures  24  and  25)  spread  widely  because  of  various 
earthquake  magnitudes  and  distances. 

47.  Damage  to  structures  in  the  epi'-entral  area  is  generally  more 
severe  on  soft  sites  than  on  hard  sites,  based  on  past  experience  and 
observations.  However,  this  study  showed  the  seismic  total  intensity 
(seismic  energy)  at  soft  sites  to  be  lower  than  at  hard  sites.  Thus, 
the  degree  of  damage  to  structures  does  not  correlate  with  the  seismic 
total  intensity  in  the  epieentral  region.  Furthermore,  the  predominant 
frequencies  at  soft  sites  are  in  the  range  of  0  to  2.5  Hz;  the  seismic 
energy  in  this  low-frequency  range  deserves  further  investigation  as  it 
relates  to  structural  damage. 
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PART  VI:  SUMMARY,  CONCLUSIONS,  AND  KKCOMMENDAT I ONS 


Summary 

48.  This  study  presents  the  results  of  a  statistical  analysis  of 
the  spectral  shapes  of  PSD  functions  of  0  to  10  Hz  for  421  ground  accel¬ 
erograms  from  89  earthquakes,  mostly  in  the  western  United  States  and 
Japan.  The  421  horizontal  accelerograms  recorded  on  ground  surface  level 
have  been  divided  into  groups  representing  (a)  rock  sites  (56  records), 

(b)  stiff  soil  sites  (131  records),  (c)  deep  cohesiouless  soil  sites 
(120  records),  and  (d)  soft  to  medium  lays  and  sands  (114  records). 

The  significant  earthquake  information  (earthquake  name,  recording  sta¬ 
tion,  date,  distance,  magnitude,  MMI,  and  peak  acceleration),  base  average 
a 

power  (area  under  the  PS!)  curve  lor  the  extended  record  of 

163.82  sec),  base  average  acceleration  t  ,  conversion  factor  (ratio  of 

163.82  sec  to  the  duration  of  the  selected  record,  or  strong-motion 

duration),  raw  average  power  A“  (spectral  intensity,  area  under  the 

RPSD  curve  of  the  actual  duration  or  the  selected  strong-motion  duration), 

2 

raw  average  acceleration  A  ,  corrected  average  power  A_  (about 

2  0  s 

12.5  percent  less  than  A“)  ,  corrected  average  acceleration  Ag  ,  and 

total  intensity  L  for  each  record  are  listed  in  Tables  1-4.  Values 
2  ° 

of  A  ,  which  were  not  directlv  estimated  Iron  original  (actual)  dura- 

°  2  2 
tion,  were  converted  from  A  .  Generally,  A  was  about  12.5  percent 

2  ° 

higher  than  A 

s 

49.  All  421  accelerograms  were  extended  to  163.82  sec  at  an  equal 
time  interval  At  of  0.02  sec  by  adding  a  string  of  zero  accelerations. 
Then,  the  statistical  mean,  standard  deviation,  etc.,  were  calculated 
for  the  extended  accelerogram  and  adjusted  to  a  zero  mean.  Next,  the 
PSD  and  the  area  under  the  PSD  curve  of  the  extended,  zero-mean  accel¬ 
erogram  (up  to  a  frequency  of  JO  Hz)  were  estimated.  Finally,  the  PSD 
curves  were  normalized  to  a  unit  area  to  contain  NPSD  curves.  The 
statistical  mean  and  mean  plus  one  standard  deviation  of  the  NPSD  curves 
for  the  four  site  conditions  were  established. 

50.  From  the  comparison  of  the  lour  site-dependent  NPSD  spectra. 


it  is  clear  that  there  are  differences  in  PSD  spectral  shapes  depending 
on  the  site  conditions.  Two  major  groups  are  formed  in  the  mean  NPSD 
spectra:  soft  to  medium  clays  and  sands  and  deep  cohesionless  soil 
sites  are  similar  and  form  one  soft  group;  stiff  soil  and  rock  sites 
form  one  hard  group.  The  frequency  of  2.5  Hz  (0.4-sec  period)  forms  a 
dividing  line  on  the  frequency  axis;  in  the  frequency  range  lower  than 
2.5  Hz,  spectral  amplifications  for  the  soft  sites  are  much  higher  than 
the  hard  sites;  and  in  the  frequency  range  higher  than  2.5  Hz,  spectral 
amplifications  for  the  hard  sites  are  higher.  However,  in  the  case  of 
mean  plus  one  standard  deviation  NPSD  spectra,  the  peak  amplitudes  in 
the  low  frequency  range  (<2.3  Hz)  decrease  in  the  order  of  soft  soil, 
deep  cohesionless  soil,  stiff  soil,  and  rock  sites  in  accordance  with 
the  degree  of  hardness.  This  order  seems  in  correlation  with  the  damage. 

51.  In  the  soft  group,  the  peak  amplitude  at  about  1  Hz  for  the 
soft  to  medium  clays  and  sands  is  about  14.3  percent  higher  than  the 
deep  cohesionless  soils.  Both  of  them  are  monotonically  attenuated  from 
the  sharp  peak  at  1  to  10  Hz  (except  one  hump  at  2.75  Hz  for  the  deep 
cohesioniess  soils).  In  the  hard  group,  the  largest  peak  amplitudes  for 
the  rock  sites  and  stiff  soil  sites  are  at  2.75  Hz  (0.36  sec)  and  0.8  Hz 
(1.25  sec),  respectively.  Generally  speaking,  the  energy  content  is 
spread  widely  over  the  frequency  range  of  0  to  10  Hz.  It  is  possible 
that  there  is  a  connection  between  the  largest  peak  amplitude  at  2.75  Hz 
for  the  hard  group  and  the  hump  at  2.75  Hz  for  the  deep  cohesionless 
soils  or  the  soft  group.  In  other  words,  it  could  be  said  that  2.75  Hz 
is  a  common  frequency  of  bedrock  under  the  deep  cohesionless  soils. 

52.  A  qualitative  comparison  was  made  of  the  spectral  shapes  of 
PSD  in  this  study  with  the  Acceleration  Response  Spectra  (ARS)  of  Seed 
and  Idriss  (1971),  Seed,  Ugas,  and  Lysmer  (1976),  and  Kiremidjian  and 
Shah  (1978).  There  is  general  agreement  in  the  spectral  shapes  of  both 
methods  except  those  for  rock  sites  that  the  amplitude  of  the  ARS  of 
Seed,  Ugas,  and  Lysmer  (1976)  was  lower  than  for  the  PSD  in  the  high 
frequency  range.  This  discrepancy  in  amplitudes  of  spectra  between  the 
PSD  and  the  ARS  was  due  principally  to  the  smaller  number  of  records  for 
the  latter.  The  number  of  rock  site  records  used  was  28  for  the  ARS  and 


52 


5t>  lor  the  PSD.  Thus,  the  avoiage  spectral  amplitude  value  for  the  former 
is  less  reliable  than  the  latter.  Another  dil h  rente  is  that  the  spec¬ 
tral  shapes  of  the  PSD  function  show  more  peaks  than  do  the  ARS. 

> 

53.  The  base  average  power  V  (.area)  under  the  PSD  curve,  which 

is  uniformly  distributed  on  163.82  see  of  the  extended  record,  has  been 
established  for  each  of  the  421  records.  The  raw  average  power  , 

or  any  average  power  of  strong-motion  duration  (see  Table  6),  which 
is  inversely  proportional  to  the  duration,  can  be  easily  calculated  as 
A“  times  the  ratio  of  163.82  sec  to  the  duration  of  the  selected  record. 

54.  The  approximate  linear  relationships  of  the  maximum  ground 

accelerations  a  versus  the  base  average  accelerations  '  and  the 

max 

average  accelerations  \  provide  a  set  of  sealing  curves  for  the  four 

°  ■, 

site  groups.  Since  the  is  the  PSD  spectral  intensity,  or  the 

scaling  factor,  a  PSD  spectrum  coaid  be  generated  for  any  desip, n  earth¬ 
quake  based  on  this  set  of  scaling  curves  and  the  four  standard  mean 
and  mean  plus  one  standard  deviation  MPSD  spectra. 

55.  In  tlie  final  analyses,  close  relationships  do  exist  among  Lhe 
three  parameters,  duration,  average  acceleration,  and  peak  ground  ace le ra¬ 
tion.  Duration  is  inversely  proportional  to  Lhe  square  of  the  average 
acceleration  (average  power).  The  latter  has  an  approximate  linear 
relation  with  the  peak  ground  acceleration.  It  is  apparent  that  duration 
has  a  large  effect  on  the  average  aeeeleration  or  PSD  spectral  intensity 

,  I.o.,  an  engineering  intensity  sc. tie.  It  is  also  the  scaling 
factor  for  the  normalized  standard  PSD  spectrum.  However,  the  average 
acceleration  is  a  relative  value  that  varies  witli  duration;  it  is  a 
scaling  factor  of  the  NPSD  curve.  it  might  be  useful  to  take  the  Arias 
Intensity  i  ^  as  a  standard  earthquake  intensity  scale.  The  Arias 
value  includes  the  total  duration  and  the  average  acceleration. 

('one  I  us  i  ons 

56.  The  statistical  analysis  of  421  accelerograms  shows  clear 
differences  in  spectral  shapes  for  dilferent  soil  and  geological  condi¬ 
tions.  Within  the  high-frequency  range  of  2.5  (o  10  Hz,  t  lie  spectrum 


for  the  rock  sites  contains  the  highest  energy  or  intensity;  the  spectrum 
of  the  stiff  soil  sites  is  slightly  lower  than  for  the  rock  sites;  and 
the  spectra  of  the  soft  clay  and  sand  sites  and  the  deep  co'  esionless 
soil  sites  are  lower  still  and  almost  the  same.  However,  in  the  low- 
frequency  range  of  0  to  2.5  Hz,  the  reverse  exists:  the  soft  sites  indi¬ 
cate  the  highest  energy,  the  deep  cohesionless  soil  sites  are  next,  the 
stiff  soil  sites  are  third,  and  finally,  the  rock  sites.  Generally,  the 
spectra  of  rock  sites  and  stiff  soil  sites  of  similar  characteristics 
can  be  classified  together  as  hard  sites;  the  other  two  site  types  can 
be  classified  together  as  soft  sites. 

57.  The  site  dependence  of  NPSD  spectra  have  been  established  by 

statistical  analysis  as  expected.  The  most  significant  finding  of  the 

study  is  the  approximate  linear  correlation  of  the  PGA  (a  )  and  aver- 

2  max 

age  acceleration  (X  ).  Since  X  is  the  area  under  the  PSD  curve, 

2  °  ° 

therefore  X^  can  be  used  as  a  scaling  factor  for  NPSD  spectra.  If 

a  is  given,  X  can  be  found  from  the  correlation  curves  of  a 

max  o  2  max 

and  Xq  .  The  standard  NPSD  spectrum  can  be  amplified  by  Xq  to  become 
a  design  PSD  spectrum. 

58.  The  comparison  of  the  attenuation  curves  of  the  PGA  and  the 
average  acceleration  versus  distance  of  the  San  Fernando  earthquake  of 
9  February  1971  showed  that  the  attenuation  rate  of  the  average  accel¬ 
eration  is  less  than  the  PGA  and  approximately  linear  on  a  log  scale. 

Recommendations 


59.  Further  developments  in  the  following  three  areas  are  needed: 

ji.  Generation  of  accelerograms  based  on  the  design  PSD 
spectrum. 

b.  Relationships  between  the  PSD  spectrum  and  the  response 
spectrum. 

c.  A  new  earthquake  engineering  intensity  scale  based  on 
average  acceleration  or  average  power  of  strong-motion 
duration. 
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(Cont  In.  nil) 
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(Japan) 
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Table  5 

Comparison  of  Average  Powers  Estimated  Frog  Uncorrected,  Corrected. 
and  Extended  163. 82-sec  Duration  Records 


Row 


Records  and  Duration 


*  From  Ravara  (1965). 


El  Centro  1940  El  Centro  1934 
N-S  E-W  N-S  E-W 


Taft  1952 


N-S 


E-W 


Olympia  1949 
N-S  E-W 


(1)  Ravara*  uncorrected 
Error  (%)  relative  to  (2) 

(2)  CIT  corrected  (standard 

value) 

(3)  Extended  163.82  sec 
Error  (%)  relative  to  (2) 

(4)  Final  correction 

(3)  *  0.875 

Accuracy  of  this  study 
(2)  relative  to  (2) 


30  sec 

25  sec 

25  sec 

30  sec 

30  sec 

30  sec 

Average  Power,  cm^sec-*1 

3820 

2690 

1900 

2305 

1360 

1775 

2930 

2190 

19.6 

17.6 

57.3 

65.9 

39.0 

65.9 

38.5 

59.3 

3193 

2288 

1208 

1389 

978 

1070 

2116 

1375 

3621 

2637 

1309 

1502 

1124 

1219 

2305 

1525 

13.4 

15.3 

8.4 

8.1 

14.9 

13.9 

8.9 

10.9 

3168 

2307 

1145 

1314 

983 

1067 

2017 

1334 

0.8 

0.8 

5.2 

5.4 

0.51 

0.3 

4.7 

3.0 

Example:  Conversion  Factor  •  rPcor{!  length,  Cgc-9) — 

Actual  or  selected  length  (sec) 


163.82 

30.0 

5.46  for  El  Centro,  1940 


Average  power  for  N-S  com-  _  Base* average  power  for  163.82  sec 
ponent  at  duration  30  sec  duration  (Column  10)  x  conversion 

factor 


663.2  (from  Table  2)  *  5.46 


2  -4 

3621  cm  sec  (in  3rd  column  of  Row  3) 


Final  corrected  average 

power  due  to  adding  zeros  0  0.875 


This  is  the  record  number  in  this  study;  i.e.,  Tables  1-4.  ^ 

The  values  of  strong-motion  duration  were  calculated  from  I  =  S  o  by  Vanmarcke  and  Lai  (1977),  where 

o  o  o  - 

This  was  calculated  from  the  values  of  I  of  Vanmarcke  and  Lai  (1977). 


Statistical  Characteristics  of  Earthquake  Ground  Motion 
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Table  8 

Correlation  of  rms  Intensity  Versus  MM1 

Upper  bound  of 
of  the  sum  of  2 

2 

the  site-dependent  rms  intensity,  cm/sec 
-horizontal  variances) 

(square  root 

Site  Conditions 

[MI 

Hard  Sites 

Soft  Sites 

:n 

400-550 

250-400 

:i 

285-395 

170-280 

r 

205-275 

120-200 

X 

145-190 

84-140 

'III 

105-140 

59-100 

In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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Errata  Sheet 
No.  1 

SITE  EFFECTS  ON  POWER  SPECTRAL  DENSITIES 
AND 

SCALING  FACTORS 
Miscellaneous  Paper  GL-81-2 
July  1981 

1.  Page  29,  Figure  9:  Change  the  horizontal  scale  of  0.1,  1,  and  10  to 
1,  10,  and  100. 

2.  Page  31,  Figure  11:  Change  the  horizontal  scale  of  1,  10,  and  100  to 
0.1,  1,  and  10. 

3.  Page  55,  References:  Line  16,  change  Keon  County  to  Kern  County. 
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